基于GPS和水准数据的微粒群算法联合反演祁连山北缘断裂三维滑动速率     

刘杰 《大地测量与地球动力学》2019,39(9):906-909

将微粒群算法与位错理论模型相结合，采用中国地壳运动观测网络提供的青藏高原地区2001~2004年GPS测量数据和2000~2006年水准测量数据，通过常规定权和附有相对权比的方法对祁连山北缘断裂的三维滑动速率进行联合反演，并与蚁群算法反演结果进行对比。结果表明，微粒群算法收敛速度快、稳定性高，结合经典位错理论模型，是一种可以有效求解断层三维滑动速率反演问题的优化算法，在大地测量反演领域极具应用潜力。  相似文献   

Gully erosion spatial modelling: Role of machine learning algorithms in selection of the best controlling factors and modelling process     

《地学前缘(英文版)》2020,11(6):2207-2219

This investigation assessed the efficacy of 10 widely used machine learning algorithms (MLA) comprising the least absolute shrinkage and selection operator (LASSO), generalized linear model (GLM), stepwise generalized linear model (SGLM), elastic net (ENET), partial least square (PLS), ridge regression, support vector machine (SVM), classification and regression trees (CART), bagged CART, and random forest (RF) for gully erosion susceptibility mapping (GESM) in Iran. The location of 462 previously existing gully erosion sites were mapped through widespread field investigations, of which 70% (323) and 30% (139) of observations were arbitrarily divided for algorithm calibration and validation. Twelve controlling factors for gully erosion, namely, soil texture, annual mean rainfall, digital elevation model (DEM), drainage density, slope, lithology, topographic wetness index (TWI), distance from rivers, aspect, distance from roads, plan curvature, and profile curvature were ranked in terms of their importance using each MLA. The MLA were compared using a training dataset for gully erosion and statistical measures such as RMSE (root mean square error), MAE (mean absolute error), and R-squared. Based on the comparisons among MLA, the RF algorithm exhibited the minimum RMSE and MAE and the maximum value of R-squared, and was therefore selected as the best model. The variable importance evaluation using the RF model revealed that distance from rivers had the highest significance in influencing the occurrence of gully erosion whereas plan curvature had the least importance. According to the GESM generated using RF, most of the study area is predicted to have a low (53.72%) or moderate (29.65%) susceptibility to gully erosion, whereas only a small area is identified to have a high (12.56%) or very high (4.07%) susceptibility. The outcome generated by RF model is validated using the ROC (Receiver Operating Characteristics) curve approach, which returned an area under the curve (AUC) of 0.985, proving the excellent forecasting ability of the model. The GESM prepared using the RF algorithm can aid decision-makers in targeting remedial actions for minimizing the damage caused by gully erosion.  相似文献   

基于最大球算法定量表征土体空间孔隙网络     

解朝阳  张巍  姚东方  王坚  朱雨辰 《工程地质学报》2020,28(1):60-68

孔隙网络控制着土体渗流、排水固结与基质吸力等重要工程性质。本文介绍了多孔介质孔隙网络最大球建模基本原理与算法。以显微CT扫描振捣干法生成的南京粉砂试样为例,采用最大球算法建立了试样三维重构模型表征单元体(REV)的空间孔隙网络球棍模型,计算得到了样品REV尺度的孔隙网络参数,统计发现,孔隙半径、喉道半径、孔隙配位数、孔隙截面形状因子、喉道截面形状因子与喉道长度等孔隙参数均近似服从正态分布,孔隙体积近似服从衰减型指数分布。孔隙与喉道半径分别分布在100μm与65μm以内,两者数学期望分别为40. 0μm与18. 0μm;配位数分布在25以内,数学期望为5. 1;孔隙与喉道截面形状因子分别分布在0. 01~0. 04与0. 01~0. 05的区间内,两者数学期望分别为0. 019与0. 033;喉道长度分布在100~800μm以内,数学期望为292. 22μm。同时发现,样品中体积小于1. 5×10^7μm^3的小孔隙数量超过90%。本方法可应用于土体细观孔隙结构的定量表征。  相似文献   

基于CUDA的地震相干体并行算法   总被引：1，自引：0，他引：1  

张全 《地质与勘探》2020,56(1):147-153

相干体技术在地震勘探资料解释方面得到了广泛的应用,由于相干体技术处理的对象是三维地震数据体,所以算法运算时间较长。为了缩短解释周期,本文充分发挥GPU并行计算优势,对C3相干体算法进行并行化分析。从硬盘读取数据到GPU上计算相干值并写入硬盘的整个过程进行分析,剔除了冗余数据的读取,完成了C3相干体算法的并行化设计与实现。最后分别对串行算法与并行算法进行性能测试,结果表明本文设计的并行算法在保证精度的前提下达到了16倍左右的加速比,对加快地震资料解释具有重要意义。  相似文献   

太原盆地及周边地区双差层析成像   总被引：2，自引：0，他引：2       下载免费PDF全文

王长在  吴建平  杨婷  王未来  范莉苹  房立华 《地球物理学报》2018,61(3):963-974

山西断陷带位于鄂尔多斯与华北地块交汇处，是我国著名的历史强震活动带之一，尤其是断陷盆地中部区域，6级以上地震频发.本研究旨在揭示忻定盆地中南部、太原盆地及临汾盆地中北部交汇处的深部结构特征，分析盆地的形成和演化，探讨该区域孕震环境.利用山西地震台网观测数据以及固定地震台站融合流动台站得到的地震数据，共7455个地震事件，采用双差层析成像方法，反演得到了太原盆地及周边地区的三维P波速度结构及精定位结果.层析成像结果显示，忻定盆地、太原盆地的中上地壳为明显的低速异常，被高速的石岭关隆起隔开，深部结构特征相对简单.太原盆地、临汾盆地及灵石隆起之间的深部结构特征较为复杂，反映了两个盆地演化过程的复杂性.穿过忻定盆地的速度剖面显示，在中地壳存在明显的低速异常体，且大部分地震都发生在该低速体上方；穿过太原盆地北部的剖面显示，该区域在20~25 km深度范围内有较密集的地震分布，并勾画出交城断裂呈犁形的断层特征；穿过太原盆地中部的剖面显示，太原盆地自西向东沉积层逐渐减薄；穿过临汾盆地的剖面揭示，汾东断裂在浅部倾角较陡，随深度增加倾角逐渐变小，倾向向东.  相似文献   

基于分区间泛化Kriging近似模型的船舶局部结构形状优化     

张干锋  王德禹 《海洋工程》2018,36(6):46-57

为提高Kriging近似模型在船舶结构性能多维度响应预测方面的适用性,对常规Kriging近似模型进行分区间泛化改进:一是对设计样本点的各个维度(分量)进行划区,并在每一个分区间内采用最优拉丁超立方(OLhd)取样;二是引入比例系数w1组合高斯与指数型相关函数提高模型对数据的泛化能力,从而在每个划分的区间内建立泛化的Kriging近似模型。通过SCH测试函数,验证了构造的分区间泛化Kriging模型有效性。结合参数化建模和改进Kriging近似模型对某过渡肘板和舱口角隅边界进行形状优化,优化第一阶段由改进的近似模型通过多岛遗传算法得到全局初步的最优解,第二阶段在初步解的基础上缩小优化变量范围,由少量的FEM计算即可在小范围内搜寻到精确的最优设计变量。结果表明:分区间泛化Kriging近似模型在预测多维度响应时较常规Kriging模型预测精度更高;分阶段的形状优化流程在保证极小误差和缩小计算成本的情况下能够得到理想的应力分布和重量优化结果,有助于船舶结构的轻量化研究。  相似文献   

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术   总被引：1，自引：0，他引：1       下载免费PDF全文

刘财  符伟  郭智奇  刘喜武  刘宇巍 《地球物理学报》2018,61(6):2589-2600

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究，针对页岩黏土含量高、层间微裂缝发育等特点，利用Backus平均理论描述页岩黏土矿物弹性参数，利用Chapman理论计算与水平微裂缝有关的VTI各向异性，并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法，并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息，并以测井纵、横波速度作为约束条件建立反演的目标函数，同时利用粒子群算法进行最优化搜索.计算结果表明，基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高，并且在一定范围内变化，这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关，同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关，而与横波各向异性参数γ相关性较小.另外，页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.  相似文献   

2017年四川九寨沟M_S7.0强震的余震重定位及主震震源机制反演          下载免费PDF全文

梁姗姗  雷建设  徐志国  徐锡伟  邹立晔  刘敬光  陈宏峰 《地球物理学报》2018,61(5):2163-2175

本文采用双差定位法对2017年8月8日至10月31日期间四川九寨沟M_S7.0主震及5200个余震序列进行相对定位，得到4036个重定位地震事件.采用中国区域地震台网观测到的宽频带垂直分向波形数据和W震相反演方法，得到了主震震源机制解.重定位结果显示，余震序列分别沿NNW和SSE两个方向扩展，展布长度约58 km，且这些余震主要集中在22 km深度之上.余震分布的另一个重要特点是具有分区特性，即在主震NNW方向约5 km处存在明显的西北和东南两区余震活动分界线；西北区的余震由深至浅具有较好连续性，而东南区却在约10 km深度处存在不连续性.余震分布的这种分区特征，说明九寨沟地震震源区的地壳结构存在强烈的不均匀性.余震分布与主震破裂特征的一致性，证实了我们定位结果的可靠性.主震的震源机制解展示出节面Ⅰ的走向/倾角/滑动角分别为246°/83.7°/-177°，而节面Ⅱ的走向/倾角/滑动角为155.7°/87.1°/-6.3°，最佳质心深度为15.5 km，矩震级M_W为6.5.根据余震分布较为垂直和主震震源机制解两节面的倾角均在80°以上，并结合野外地质调查结果，推测此次九寨沟地震为与节面Ⅱ参数相近的一次高角度的左旋走滑型事件.  相似文献   

基于F-范数的不确定性平差模型的解算方法     

鲁铁定  朱国红 《大地测量与地球动力学》2018,38(6):557-561

为提高基于F-范数的不确定性平差模型的解算效率，给出直接迭代算法进行参数估计。该算法无需SVD，解算过程简单且易于编程计算，同时给出迭代不收敛时的SVD-解方程算法。二元线性拟合及沉降观测AR模型的算例结果表明，这2种算法正确可行，与SVD-迭代算法具有等价性。当迭代收敛时，宜使用直接迭代算法，收敛速度更快，解算效率更高；当迭代不收敛时，可釆用SVD-解方程算法。  相似文献   

Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*          下载免费PDF全文

Yuyang Tan  Haijiang Zhang  Junlun Li  Chen Yin  Furong Wu 《地震科学（英文版）》2020,33(4):177-193

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.  相似文献